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Abstract 



Stirling's formula, the asymptotic expansion of n! for n large, or of r(2;) for z -^ oo, 
is derived directly from the recursion equationr(z + 1) = zT{s) and the normalization 
condition F f ^ j = ^/tt. 



Stirling's formula is an asymptotic expansion of the gamma function T{z), valid as 2; — > oo 
with I argz| < tt. Since T{z) is the analytic extension of (n — 1)!, it satisfies the recursion 
equation 

T{z + l) = zTiz). (1) 

We shall derive Stirling's formula from ([T]) and the condition 

r (i) = y?. (2) 

which makes the solution unique. 

Taking the logarithm of each side of ([1]), and writing f{z) = logr(2;), yields 

fiz + l)-f{z) = \ogz. (3) 

Replacing f{z + 1) in ([3]) by its Taylor series around z leads to 

n=l "■■ 

We integrate (jl]), introducing a constant of integration log C, and we call the right side fo{z): 

E ^ = ^log^-^ + logC^/o(^). (5) 

Keeping the undifferentiated term with n = 1 on the left, and moving the other terms to the 
right, gives 

/W = /oW - E ^^^- (6) 

Equation ([6]) is a differential equation for / of infinite order, which we solve by successive 
substitution. First we approximate / by /q. Then we use /o in place of / in the sum to get 
a second approximation to /. Then we substitute that for / in the sum, and so on. This 
will lead to an expansion for / of the form 

CO 

/W-E^tf^^)- (7) 

j=0 



We expect ([7]) to be an asymptotic expansion of f{z), valid as z -^ oo. Therefore we have 
used the sign ~ of asymptotic equahty. 

To determine the coefficients Cj in ([7]) we substitute ([7]) into IQ and equate coefficients 
of /o . This yields 

Equation (IHl) determines the Cj successively, starting with cq = 1, ci = —1/2, etc. To solve 
it in general we set Cj = Bj/j\ in ([8]). The resulting recursion equation for Bj is exactly that 
satisfied by the Bernoulli numbers [eq. (6.79), p. 270 of Graham, Knuth and Patashnik [1]]. 
Thus Bj is the j-th Bernoulli number. 

The expansion ([7]) of / involves the derivatives of /q. Since /g = log 2; it follows that 
f^^^ = {-iy{j - 2)\/z^~^ for J > 2. For j > 2 with j odd, Bj = 0. Therefore by setting 
j = 2n for j > 2, and using the definition f{z) = logr(z), we can write (JTj) as follows: 

log r(^) ~ (^ - i) log . - ^ + log c + 1 %iT^ - (9) 

Exponentiating both sides of iQ gives 

{2n-2)\B2n' 



T{z) ~ Cz"" 2e "^exp 
To find C we first prove that 



S i^ny.z 



2n-l 



(10) 



V 2/ 22"r(n + l) ^ ^ 

For n = this reduces to ([2]) , so it is true. It follows for ^ = n by induction on n using (JT]). 
When we use the leading term in flTOj) for T{n) in fITT]) . we find that 

C = V2^. (12) 

This completes the determination of the expansions iQ and fITOl) of logF and of F. With 
C = -\/27r, Q) agrees with eq. (6.1.40) of Abramowitz and Stegun [2]. When the exponential 



in (TTO!) is expanded in a power series, it agrees with eq. (6.1.37) of Abramowitz and Stegun 
^, which is Stirhng's formula. 
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